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Abstract. In this work we analyze the dynamical behavior of the collision between two clouds of fermionic 
atoms with opposite spin polarization. By means of the time-evolving block decimation (TEBD) numerical 
method, we simulate the collision of two one-dimensional clouds in a lattice. There is a symmetry in 
the collision behaviour between the attractive and repulsive interactions. We analyze the pair formation 
dynamics in the collision region, providing a quantitative analysis of the pair formation mechanism in 
terms of a simple two-site model. 



1 Introduction netic states in ultracold gases have been considered in 

the experiment Q]. Indeed, to understand the formation 

g • In recent years, ultracold gases have become an unpar- and properties of strongly correlated states of Fermions 

^ alleled tool for simulating condensed matter systems [1] it is essential to study also their dynamics. In condensed 

_ > ! and to explore the properties of paradigmatic condensed matter physics much work has been devoted to the deter- 

, matter models. Ultracold gases allow an unprecedented niination of the ground-state properties of the Hubbard 

tunability of the system parameters and dimensionality Hamiltonian (see e.g. 0), while its dynamical behavior 

For instance, the experimental reahzation of the Fermionic has been explored to a much lesser extent. Nevetheless, 

Hubbard Hamiltonian [IE] represents an important effort in the recent past, the versatility of ultracold atomic sys- 

towards studying phenomena underlying high tempera- tems has led, from a numerical and theoretical point of 

ture superconductivity and its connection with antifer- view, to approach the analysis of the dynamics in such 
romagnetism in cuprates and iron pnictides. Ferromag- 
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systems |6l7l8l9ll0lllll2ll3ll4j . leading to a revived in- 
terest in the unitary evolution of closed quantum systems 
|15| . Recently, an interesting experimental investigation of 
spin dynamics in a system of colliding Fermi gas clouds system. Spin up and down gases in a lattice 

was reported [IS], closely related to the topic of this arti- ^0^^^"^^ in two separate harmonic potentials. At time 

zero the harmonic potentials are removed and the gases 
expand and collide. Here, L and R mark the two central 

In this work, we simulate the collision of spin polarized 

sites where the expanding gases meet and which we use 

gases using the time-evolving block decimation (TEBD) 

for the two-site model. 

algorithm |17j . We are interested in the coUisional proper- 
ties of two clouds with opposite spin polarization (denoted discuss the connection between our resuhs and the resuhs 
up and down hereafter). Initially, the two polarized gases in the low-temperature strongly interacting regime of the 
are trapped by separate harmonic potentials. At t = we recent experiment [H]. 

turn off the harmonic traps, allowing the clouds to expand In section [2] we describe the system in detail. In sec- 

and colhde with each other, in complete analogy to what tion[2]we elaborate on the TEBD numerics and in section 

has been done in the experiment reported in [IS], except H explain the Hubbard Dimer two-site model. Finally, in 

for the presence of the lattice. section [5] we compare the results of the numerics and the 

n 1 • 1 1 • • -1 model. 

Both attractive and repulsive interactions between the 

species are considered. Counterintuitively, the physics of 

the collision is independent of the sign of the interaction. 2 The system 

From an intuitive point of view, one might expect that. 

We consider one-dimensional spin up and spin down gases 

in presence of attractive interaction, the particles would 

confined in separate potentials in space. In addition to 

merge in a gas of pairs, while bouncing off for repulsive in- 

the harmonic trapping potentials, there is a lattice poten- 

teraction. However, the actual quantum unitary dynamics 

tial (see Fig.[T]). The system is described by the Hubbard 

is different. Of particular interest is the pair creation dur- 

Hamiltonian with a harmonic potential: 

ing the collision for both signs of the interaction strength. 
In the present article, it is shown that the mechanism for 

pair creation in the collision is explained by a two-site ^ 4- + ^ X]i(^t *) "^t 

analysis previously introduced by us [T5]. It is also shown +^4-Si(^4- *) 'JJ2ia=^,i'^i(7^i-+icr + h.c, 

that the qualitative dynamics of the collision fundamen- where the c\ ^ operator creates a spin down particle at 
tally arises from this Hubbard Dimer model. Finally, we lattice site i, c; + annihilates an up particle at lattice site 
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i, fii,^ = cl^Ci^i-, J is the hopping matrix element, U is Hamiltonian and the latter determines the time evolution 

the interaction strength between the spin up and down when the wavefunction is acted on by e~*^', i.e. by solving 

particles, and V and C are the spin-dependent harmonic the time-dependent Schrodinger's equation. The Schmidt 

trapping strength and position of the trap center, respec- number for TEBD (describing the numerical truncation) 

tively. in the simulations is _r = 150. 



The initial state of the simulations is the ground state 

of the Hubbard Hamiltonian ^ with V\ ^ ^ 0.02, 

„ rr, r 1 ^ 1 r ml ,1 • ^ ^ ■ i i Up aud 20 dowu particlcs, i.e. Ni — 20, — 20. The in- 
= 50.5, and Cj. = 100.5. Then, this state is evolved f : ^. ^ \ 



For the simulation parameters we choose to have 20 



tcraction is varied so that we run the simulation for inter- 
actions ^ = 0,±1,±3,±5,±7,±10,±15 and ±20, where 
the interaction is expressed in the units of the hopping 
J. Above and from now on, variables are expressed in the 
units of hopping, and J = 1 has been chosen in the nu- 
merics. Note that in our convention negative values of U 

(TEBD) algorithm [17] to model the collision in one di- , ,, ,. . , , „ ,. mN mi 

' ' represent attractive interaction (see Equation (UJ). ihe 



with the otherwise same Hamiltonian, except V(i)t = 
Vii)i = 0.0. 

3 TEBD numerics 

In this work, we use the time-evolving block decimation 



initial trapping strengths are chosen to be -j = 0.02. We 
consider the temperature T = 0. At time zero, we release 
the traps (i.e. change V from 0.02J to O.OOJ) and let the 
two clouds expand, keeping the interactions on during the 
expansion. 



mension and in a lattice. The only approximations in- 
volved in the TEBD numerics are in the Suzuki- Trotter 
expansion and in the Schmidt truncation of the Hilbert 
space [17] . TEBD is an essentially exact numerical method 
in the sense that it does not use mean-field approxima- 
tions, and the errors due to the approximations above are 

controllable. For the details of the TEBD numerics, see We run the time evolutions up to the time t = 25 j, 
[It] . as then the outer edges of the clouds have hit the edges 

The initial state for the time-evolution is calculated us- of our finite system (the lattice size L = 150), and we 
ing a ground state algorithm for TEBD [17j . We have used are not interested in the unphysical edge collision dynam- 
two different algorithms, one for the time evolution and ics. However, the maximum speed for propagation of the 
another for the ground state, and both of the algorithms distrurbance caused by the collision with the edge is 2 J 
employ TEBD. The ground state calculating algorithm is (due to the lattice dispersion). As we are interested only 
called the "imaginary time evolution algorithm" and the in what happens in the collision center, we can run the 
time evolution one is called the " real time evolution algor- simulations a little longer than the time when collision 
tithm" . The former solves the ground state of the Hubbard with the edges occurs. 
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As a result of the TEBD simulations, we obtain the 
density profiles of up and down particles, ni^^{t), ni^i{t). 
Importantly, also the density of doublons, ni^-fi{t) is ob- 
tained: 

n^uit) ^W|cItC^t4iC^i|^(^) >, (2) 

where <> denotes the quantum mechanical expectation 
value, and (l>(t) is the wavefunction. Elaborating on the 
definition of doublons, the doublons are excitations of the 
form c,|^c|j^|0) and the single (unpaired) particles are de- 
fined as c|^|0) (ct =t, I), where |0) is the state represent- 
ing an empty lattice site. The local number of doublons is 
given by ni-\-i{t), while the number of unpaired (up) par- 
ticles is given by f^""(^) — ni^{t) — ni^i{t). Now, before 
moving on to the results of the simulations let us discuss 
the theoretical model that we employ in order to explain 
the results. 

4 The Hubbard Dimer two-fluid model 

Previously, we have developed a Hubbard Dimer two-fluid 
model to explain the dynamics of expansion in a ID Fermi 
gas in a lattice [T5] . There we considered the expansion 
of an interacting two-component gas which was initially 
set into a band insulator state, related to the experiment 
of [in]. As will be shown in this article, the same model 
explains dynamics of the collision of ID polarized cases. 
Below, we will go through the derivation of the Hubbard 
Dimer model in the case relevant for this problem (it is ba- 
sically the same analysis as done in the online supporting 
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material of |18j but is included here, and done in greater 
detail, for clarity). 

We assume that the important dynamics occurs in the 
collision center, i.e. the two lattice sites in the center where 
the expanding polarized gases meet, see Fig [TJ The spin 
basis for a single lattice site is |0 > (empty lattice site), 
I t>: I i>: a-nd I ti>- Let us now assume that in the 
two-site collision center, just after the first particles have 
entered the system from the rest of the chain, we have the 
state \(I>{t — 0) >= I t!i>- We want to determine how 
this state evolves into a doublon state, | ti, > or |0, 'fl> 
as a function of time. In order to do that, we solve the 
two-site system with the Hubbard Hamiltonian exactly 
by diagonalizing it. The two site system and its solution 
is in general called the Hubbard Dimer (see e.g. |20|). 

We have a 2-particle basis, and the Hubbard Hamilto- 
nian conserves the number of particles. Due to anticom- 
mutation relations, a given order for the application of 
fermionic operators must be chosen: 

lt,;>=clt4jo>, I ;,t>= 4^4^10 >, (3) 
in,o>=ci/ijo>, |o,n>=4t4jo> . (4) 

The Hamiltonian is 

H = Hj + Hint 
Hj = -J^ 4,<tC2 a- + h.C. 

Htnt ^ n^-^Uii. (5) 

i=l,2 

Then 
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The Hamiltonian can be rewritten in a basis where it as- 
sumes a block-diagonal form 

= j(in,o>-|o,n>), (6) 



H 



bl — 












u 














-2J 





2J 


u 



(11) 



Hj\U> = Hjc\^4^\0> 

= j(io,n>-in,o>) (7) 

and, obviously, 

Hj\n,0> ^ Hjcl^cl^\Q> 

= -j (4tcit4t4 JO > +4;cu4^4 JO >) 
= j(it,;>-u,t>), (8) 



analogously, the hopping "perturbation" part of the Hamil- 
tonian assumes the form 



H 



J,bl 



























-2J 





2J 






(12) 



Hj\o,n> = Hj 4^4^10 > 

= -J (4^c2t4t44,lo > +4iC24,4t4;lo >) 
= j(it,;>-u,t>)- (9) 

Hence in the 4-dimensional Hilbert space of the N — 2 
particles Hubbard Dimer, with the choice of the basis 
given by Eq. (representation of Hint is trivial) , H has 
the following representation 

J -J 

-J J 

J -J U 

-J J U 



This representation corresponds to the following basis vec- 
tors 

\T>=^i\t,i>+\iA>) 

\D- >= -i=(in,o> -|o,n>) 
i^>-^(it,;>-u,t>) 
i^+>=^(in,o>+|o,n>). (13) 

If the lower block of H is diagonalized, one obtains the 
following expression for the eigenvalues 



A± = U/2 



1±A 1 



16J2 

C/2 



Defining 



H 



a± 



U 
4J 



1± Wl 



16J2 



(14) 



(15) 



(10) 



the eigenvectors can be written as 



v± >= 



1 



|5 > +a±\D+ >) 



1 + al 



(16) 
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Hence the full spectrum of the dimer is given by 
A_ = C//2 1 - ^|l 



16J2 



(< 0) 



Ao = 
\+ = C//2 



1 + 



C/2 



(> U) 



\v- > 
\T > 



\v+ > 



(17) 



Now, in the case of the problem in hand we have initially 
the state 

|(/.(t = 0)>-|t4>=^(l^>+|r>) (18) 

which we need to express in eigenstates of the Hamilto- 
nian. Expressing jiS > as a superposition of |w+ > and 
If- > gives: 



1 14>= > > +1^ >)' (19) 



where we have denoted 



(^1+^)0!=, 



(20) 



(a-l_ — Q!_) 

Now let us determine the number of doublons in the left 

site, n^^{t) —< (j){t)\n^^\<l){t) > given by 



<t,i \e « h^^e " I t,i> • 
Calculating this gives 



(21) 



16+i 



[1 - cos{VlP+TKPt)], (22) 



which determines the time dependence of the doublons in 
the problem, completing the analysis. 

Now, we need to return to our initial assumptions. 
Equation (gS]) holds for the initial state |<?(i = 0) >'^= | t 



,4,>. If more unpaired particles did not enter the system. 
Equation ([^^ would predict that we simply see oscilla- 



tions in the doublon density with the frequency \/U^ + 16J^ 
and amplitude — However, when the two gases col- 
lidc, during the collision more unpaired particles enter the 
system, going to the state | f, i>- We make the hypothesis 
that the number of doublons in the central sites n^^ (t) can 
be determined by the following short time approximation: 

*[1 - cos{^/U^ + 16J2(r - t'))] * sin^iJ * t) (23) 

where C denotes the two-site collision center. Equation 
(|23p contains the time evolution of pairs given by the 
Dimer problem. Equation (j22p . In addition, it takes into 
account the number of unpaired particles at the sites sur- 
rounding the two central sites, n^n{t ), which can change 
during the collision. The particles from these neighbour- 
ing sites tunnel into the central two sites as described by 
the term sin'^{J * t ). Particles tunnelled into the central 
two sites at time t start the Dimer dynamics at that time, 
thus the shift t — t in the cosine term. The factor of two 
comes from spatial symmetry. 

Finally, let us discuss the limitations of the above anal- 
ysis. Equation (|23l) we assume to hold for short times, 
since for longer times one needs to take into account un- 
paired particles as well as pairs tunnelling out of the two- 
site dimer. Incorporating these two into the analysis is 
somewhat meticulous, but as we are, for now, more inter- 
ested in whether the dynamics is fundamentally explained 
by the Hubbard Dimer model than whether the model 



Jussi Kajala et al.: Collision of one dimensional ( 

dynamics can be analytically solved at long times, we re- 
strict ourselves to the short time limit. The short time 
limit means the times when the change in n^ni^) due to 
tunnelling into the central sites and due to pair creation is 
small, and the pair tunnelling is neglibile. The change in 
n^„{t) is small when the sin'^{JH) — jg close to 

zero. That holds when t << ^j. This is shorter timescale 
than the pair tunnelling timescale for all U and therefore 
is the limiting timescale in our short time analysis. 

This completes our analysis of the Hubbard Dimer. Let 
us next use the obtained analytical results to explain the 
numerics. 

5 Results and Discussion 

The square roots of the up density profiles \/ during 
the collision of the oppositely polarized gases are shown 
in Figures [2] - [7] for different interactions. We are plotting 
the square roots of the density distributions since they 
highlight low density features which are important for the 
analysis in the case of pairs. Note that due to symmetry 
the density profiles of down particles n-f,i{t) are mirror 
images of \/nfJ(t) with respect to the collision center. 

Looking at Figures [5] - [5] one interestingly observes at 
every interaction \U\ that there is a C/ o —U symmetry 
in the collision. The symmetry holds for all observables 
we determined. In the case of both the attractive and 
repulsive interactions the clouds bounce back from each 
other. One way of explaining, for large \U\, this somewhat 
surprising behaviour is noting that the lattice dispersion 
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Lattice site 

Fig. 2. The square root of the density profile of up par- 



ticles, y/nf~i(F), the interaction U = 0.0. y/ni~i(J) is not 



shown as it is symmetric to -y/ n-f^i{t) with respect to the 
central lattice sites 76-77. (See online for colour). 
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Fig. 3. The same as figured! but U = —1.0. 
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Fig. 4. The same as figure [H but U = +1.0. 




Lattice site 

Fig. 5. The same as figure [21 but U — —5.0. 
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Fig. 6. The same as figured but U = +5.0. 




Fig. 7. The same as figured but (7 = — 15.0. 
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Fig. 8. The same as figured but U = +15.0. 

limits possible kinetic energy in the single band Hubbard 
model, the maximum energy being 4J. When \U\ is large, 
the large energy mismatch between a paired state and 
a non-paired one suppresses the probability of creating 
a pair from the colliding initially unpaired clouds, as re- 



flected in the Lorentzian form of the amplitude 



-r- ui 



2 1+-^., 

' 1B.J-' 

Equation (E^ . Therefore the polarized clouds are reflected 
from each other. Below, we shall compare the amount of 
doublons created in numerics to the predictions of the 
Hubbbard Dimer. Noting the U ^ ~U symmetry we will 
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Fig. 9. The square root of the density profile of dou- 



blons, ^Jn-<f^i{t), the interaction is \U\ = 1.0. (See online 
for colour). 
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Fig. 10. The same as figured but \U\ = 5.0. 

henceforth in the discussion denote interactions with ab- 
solute values. 

In Figures [9] - [TT] we plot the square root of the den- 
sity of the doublons nfi,i{t) for interactions \U\ — 1.0, 
\U\ ^ 3.0, and \U\ = 5.0. Looking at Figures d - El we 
see that doublons are indeed initially created in the colli- 
sion center, and then they spread and possibly dissociate 
back to unpaired particles. Intriguingly, we see oscillations 
at the collision center. To examine the behaviour of the 
doublons better, let us plot the total number of doublons 
n'^^*"-'' (t) , given by 



-.Total 



(24) 
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g. 11. The same as figure [HI but 
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Fig. 12. n^f^t) as obtained from TEBD numerics. 
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Fig. 13. n9.{t) as obtained from TEBD numerics. 



where L is the lattice size L = 150. We also plot the total 



number of pairs in the collision center sites n£, (i), given 



by 



where the sites L and R are the two dimer sites what we 
call "the collision center", see Figure [TJ The quantities 
n'^^*°''-{t) and n^^{t) are shown in Figures [T^ and [T^ for 
high interactions. 

Figures [12] and [13] tell about the dynamics of pair cre- 
ation both at the collision center and in total. For short 



0.02 



0.015 



o<- 0.01 



0.005 



-TEBD 

-Short time Hubbard Dimer 




Time (1 / J) 



Fig. 14. Initial growth of n^'^{t) as predicted by the 
Hubbard Dimer model compared to the TEBD numerics, 
\U\ = 5.0. 
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-TEBD 
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Fig. 15. Initial growth of n^^{t) as predicted by the 
Hubbard Dimer model, compared to the TEBD numer- 
ics, \U\ = 10.0. 

times just after the collision, the results seen in Figure [T51 
should be compared to Hubbard Dimer predictions, c.f. 
Equation (1^^ . These are shown in Figures [TJ] and [TS] for 
iji = 5.0 and ^ = 10.0. 

The plots in Figures[T^and[T51have been obtained from 
Equation (1231) . assuming the density time-dependence to 
be the shape of a square pulse, i.e. n^ni^) = 0.2. In the 
numerics, the incoming cloud shape is not square but more 
like a Gaussian with the height 0.3 and the half-width at 
half maximum of 10 lattice sites. However, assuming a 
square shape simplifies the analysis significantly and is a 
reasonable approximation, since the high density part of 
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the cloud is anyway the one that contributes the most to and down particles. When the collision progresses in time, 

the pair creation. more and more sites further away from the collision site 

Elaborating on Figures [H and [HI in these Figures is ^i^l ^ave have both up and down particles. Let us now 

shown the predicted growth of n^^{t) from the beginning define reaction center R as all the lattice sites which have 

of the collision. The differences between the predictions nonzero population of both up and down particles. More- 

and the TEBD date are likely to be due to the assump- o^^r, we define the reaction edge sites Edge{t) to mean 

tion of a step function for the incoming unpaired density ^^'^ l^st sites that have both up and down particles, when 

pulse. As said above, the actual incoming polarized cloud counting from the two sites of the collision center. These 

shape is more like a Gaussian, and thus we have neglected edge sites change as a function of time as the collision 

the small density increase in the beginning of the colli- progresses. 

sion. Indeed, initially the step function underestimates the Now, unpaired particles will tunnel into the reaction 

Gaussian, and at later times (when the short time approx- center from the sites which are adjacent to the edge sites, 

imation starts to break down) it overestimates it. The density of unpaired particles in these sites is denoted 

Looking at Figure M one sees the beginning of the analogously to the short time analysis. However, 

pair formation when the two clouds collide, at t ^ 7j. ^^^^ ^^^S*'^ considered, we need to take into ac- 

rrii • • ■ J -J- ■ 1 J ri i count also unpaired particles tunnelling out of the reaction 

Ihe maximum m the pair density is reached after r w ft- a 

o 1 r 4. i -c 11- ■ -\/t 1 center. They tunnel out from the edge sites, at which the 

3 J from the start oi the coUision. Moreover, we see large ° ' 

1 -n • ,1 • J -ii • J oi density is n^lf^'^(f). Finally, in determining the total den- 
scale osciUations m the pair density with a period w 34. un \ j & 



„ ■ tr T,!, J T-v- u J • /HToii sity of unpaired particles in the reaction center, we must 

Mow, since the Hubbard Dimer - based approximation (j23D j f t- i 

11 • rj ir j-^^TTi consider that the unpaired particles will convert into dou- 

we considered above is valid only tor times t « we ^ ^ 

, , . ,1 , , -iw • blons via Dimer dynamics. Summing up these contribu- 

cannot explain the large scale oscillations or the maximum o f 

f r+j -+1,+ + J j+ J tions, we obtain for the total density of unpaired particles 

ot amplitude using that model. Instead, we need to expand ' j f t- 

,1 1 , , • -J • i 1 rnu- ■ in the reaction center 

the sliort time considerations into longer times. Ihis is, as 

mentioned above, somewhat meticulous, but we hope to 

gain insight by formulating the general time-dependent ^uni^) = 2* {sin{Jt))'^n^^{t) 

(26) 

equations although solving them might be difficult. ~2 * {sin{Jt))'^n^^^'^{t) — n^^{t), 

In the case of longer times, we must take into account where n^'^{t) is the total number of doublons in the re- 

the fact that the dimer dynamics occurs in several other action center. Restating, the first term accounts for un- 

sites, not just at the collision center. The dimer dynamics paired particles entering the reaction center, the second 

will occur at every site which has population of both up term accounts for unpaired particles leaving the reaction 
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center and the last term accounts for unpaired particles Equations [57] - ^\ could be solved self-consistently to 

converted into pairs. Next we consider the doublons. With obtain the full time evolution predicted by the model, but 

the definitions above, we hypothesise that the growth (G) this is beyond the scope of this article. Let us instead 

and decay (D) of are given by: see if one learns something from the equations without 



^(^) - It=0 It'=t 16+^ 

*(1 - cos{VU^ + 16J2(t - t')) * n^„{t)dTdt\ 



(27) 



solving them. In the high-interaction limit, we note that 
the cosine oscillations occur at such a high frequency that 
they average out. Thus, substituting G{t) and D{t) into 
Equation [511 and evaluating the r integral: 



*(1 - cos{VU^ + 16 J2(t - t')) * flf^{t)dTdt , 
and the total number of doublons in the reaction center 



(28) n^,{t) = t.-^. nL{t)dt. (30) 



IS 



In the high interaction limit \U\ > 3, n!^^{t) in Eauation[26l 
is neglibile, as the Hubbard Dimer prefactor — ^775- makes 

" 16+% 



hf^it) = Git) - D{t). (29) 



the density or doublons produced much less than the num- 
ber of unpaired particles entering the reaction center from 
the expanding polarized clouds. Thus, in the high interac- 
It is noted that as we have defined the reaction center tion limit we obtain: 
consisting of all the sites which have nonzero population 

of up and down particles, the quantity n5 (t) is equal to ~r / ,\ _ r * =* + s 

n^°*"'(t) (unlike n^„(t) which is not equal to nlZ^^^t) . , , . , Edneu'\\^-^' 

*{2 * {sm{Jt )yK„{t ) - 2 * {sm{Jt )YK^3'^{t ))dt . 

since there exist unpaired particles outside the reaction (31) 

center). Therefore, n^\^{t) can be directly compared to the Which is our final result. The density waves of expand- 

total doublons densities seen in Figure [12] as a function of ing gases have the shape of a Gaussian, i.e. n^nit) and 

time (not to the two-site collision center doublon densities n^n^'^it) are Gaussians. To be more exact, the Gaussians 

seen in Figure [T51 that were considered above). The short have the mean 8-j (we start to coimt the time here at the 

time approximation we did initially involved 1) approxi- beginning of the collision) and half width at half maxi- 

mating the nj*^ term being neglibile in Equations (j26p and mum 10 lattice sites (these, again, depend on the shapes 

(|55|) 2) neglecting the tunnelling away of unpaired parti- of the polarized cloud). Therefore, based on the form of the 

cles, i.e. the second term in Equation ([55]) 3) considering equation [3T] we would expect to see in n^^'^°'''{t) 1) propor- 

that the relevant dynamics occurs at the two central cites, tionality to — as a function of interaction 2) Linear 

i.e. Edge{t) — L, R. increase of amplitude superimposed on a wide Gaussian - 
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Fig. 16. The height of the first peak of n^°^°-^{t) seen in 
TEBD numerics fitted to a * 



shaped increase until 8-j from the beginning of the colli- 
sion after which wide Gaussian - shaped decrease. 3) Os- 
cillations with period tt, because the only oscillating term 
in the equation has the form (sm( Ji))^ = i~cos(2jt) ^ rpj^jg 
is promising since in the TEBD data in Figure [T^] one 
sees oscillations at a period « 3j superimposed with a 
Gaussian shape increase until t = 6j + 8j = 14-j. To 
examine whether the prediction (Equation (I3ip ) for the 
interaction dependence of the amplitude of these oscilla- 
tions matches the TEBD numerics quantitaively, we plot 
the densities at the first visible oscillation peaks in Figure 

[niat i = - lOi and fit the resuh to — The result 
J J 16+^ 

of the fit is shown in Figure [1^] 

Intriguingly, we see that the fit in [1^] is very good. 
The fitting parameter a was determined to be a « 4.0 
which means that, in the long time limit, Hubbard Dimer 
dynamics indeed occur in several lattice sites close to the 
collision center. 

6 Conclusions 

We have simulated the collision of two polarized gases in 
ID and in a lattice using TEBD numerics. We found that 



there is — [/ ^ C/ symmetry in the collision and the gases 
bounce back from each other for interactions \U\ > 0.5. 
Indeed, in ID particles cannot pass each other without 
interacting. We propose that our analysis based on the 
Hubbard Dimer [18| . which is a two site model, explains 
dynamically how pairs are created and dissociated during 
the collision. This simple model explains why the gases 
bounce back: the Hubbard Dimer dynamics constrains 
the number of pairs that can be created in the time that 
the unpaired up and down particles are in contact during 
the collision. Indeed, we compared the short time Hub- 
bard Dimer analytical results to the numerics and found 
a good correspondence. In addition, we formulated time- 
development for long times (self-consistent equations) and 
were able to identify the prominent features of long time 
pairing dynamics by examining the form of the equations. 

The derivation of the Hubbard Dimer dynamics does 
not importantyly include dimensionality dependence. It is 
possible that the same dynamics works for higher dimen- 
sions, as is suggested by our earlier work |18| matching 
partly the results of a 2D experiment [19] . Interestingly, 
the simulations and analysis presented in this article could 
be mapped to the continuum case in the low density limit, 
possibly relevant for the experiment 16 . in the quantum 
unitary evolution regime; in [16 the dynamics of the col- 
lision was explained using semi-classical Boltzmann equa- 
tions, which do not describe the low temperature regime. 
However, the Hubbard Dimer dynamics in the continuum 
limit is a subject of a further study. Nonetheless, one can 
note the bouncing back - behaviour of the clouds, due 
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to interaction, both in the experiment 16 and our sim- in a fermi gas of ultracold atoms. Science, 325(5947): 1521-4, 

ulations. Moreover, our results could give a quantitative Sep 2009. 

prediction for the spin diffusivity at T = 0. It should also 5. Fabian H. L. Essler, Holger Frahm, Frank Gohmann, 

be feasible to prepare Fermi gases in ID lattice and in Andreas Kliimper, and Vladimir E. Korepin. The One- 

such sysytems our predictions could be directly tested. D^mens^onal Hubbard Model. Cambridge University Press, 

Sep 2005. 
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